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£NJ , We discuss uncertainties in the extraction of parton distributions from 

global analyses of DIS and related data. We present conservative sets 
of partons, at both NLO and NNLO, which are stable to x, Q 2 , W 2 cuts 
| on the data. We give the corresponding values of cts{M%) and the cross 

C*") ■ sections for W production at the Tevatron. 

O; 

\ The parton distributions of the proton, which are currently determined from 

O j' global analyses of a wide range of DIS and related hard scattering data, are 

^ ■ subject to many sources of uncertainty. There are uncertainties due to (i) the 

D , experimental errors on the data that are fitted in the global analysis, (ii) the 

' choice of data cuts (W cut , x cut , Qcut)' defined such that data with values of W, 

x or Q 2 below the cut are excluded from the global fit, (iii) the truncation 
, of the DGLAP perturbation expansion, (iv) specific theoretical effects, such 

as In 1/x, ln(l — x), absorptive and higher-twist corrections, and (v) input 
assumptions, such as isospin invariance, the choice of parameterization, heavy 
target corrections and the form of the strange quark sea. 

So far, attention has been focussed on the uncertainties arising from the 
experimental errors; see Refs. [1, 2] for estimates based on global NLO analyses. 
However, Ref. [3] concentrates on the remaining uncertainties, (ii)-(v); here we 
present some results from this forthcoming paper. In principle, if the DGLAP 
formalism is valid and the various data sets are compatible, then changing the 
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data that are included in the global analysis should not move the predictions 
outside the error bands. In practice this is not the case. Consider, for instance, 
the effect of different choices of x cut on the data that are fitted. Table 1 shows 
the values of x 2 f° r NLO global analyses performed for different values of a; cut , 
together with the number of data points fitted. Each column represents the % 2 
values corresponding to a fit performed with a different choice of the cut in x. 



3-cut • 





0.0002 


0.001 


0.0025 


0.005 


0.01 


# data points 


2097 


2050 


1961 


1898 


1826 


1762 


X 2 (x > 0) 


2267 












X 2 (x > 0.0002) 


2212 


2203 










X 2 {x > 0.001) 


2134 


2128 


2119 








X 2 {x > 0.0025) 


2069 


2064 


2055 


2040 






X 2 {x > 0.005) 


2024 


2019 


2012 


1993 


1973 




X 2 (x > 0.01) 


1965 


1961 


1953 


1934 


1917 


1916 


Ai +1 




0.19 0.10 


0.24 0.28 


0.02 



Table 1: The measure of stability, A- +1 , to changing the choice of x cu t- 



To obtain a measure of the stability of the analysis to changes in the choice 
of x cut , we compare fits in adjacent columns, that is with (x cu t)i+i and (x C ut)i- 
In particular, it is informative to compare the contributions to their respective 
X 2 values from the subset of data with x > (x cu t)i+i- If stability were achieved, 
then we would expect the difference A% 2 between these two x 2 contributions 
to be very small. We stress that these two x 2 contributions describe the quality 
of the two fits to the same subset of the data. Thus a measure, A* +1 , of the 
stability of the analysis is A% 2 divided by the number of data points omitted 
when going from the fit with (x C ut)i to the fit with (x cu t)i+i- For example, if 
we raise the x cu t from 0.001 to 0.0025 then Ax 2 = 2055 - 2040 for the data 
with x > 0.0025, and the number of data points omitted is 1961 — 1898 = 63. 
Thus the measure A°;°°f = 15/ 63 = 0.24, as shown in the last row of Table 1 . 

Inspection of the values of A* +1 shows a significant improvement in the 
quality of the fit each time x cut is raised until the final step when x cut is 
increased from 0.005 to 0.01, when we see that there is no further improvement 
at all. In fact, raising x cut from 0.01 to 0.02 confirms this stability. Hence we 
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conclude that x ~ 0.005 is a safe choice of a; cu t- 

After similar studies of the effect of varying Q 2 ut an d Wcut > anc * combina- 
tions of these cuts, we find a restricted domain (x cut > 0.005, Q 2 ut > 10 GeV 2 , 
VF 2 > 15 GeV 2 ) in which the obtained NLO parton distributions are stable to 
further cuts. We denote this conservative parton set by MRST03(cons). Their 
behaviour relative to MRST02 [2] is shown in Fig. 1. 
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Figure 1: The conservative partons compared to MRST2002 [2]. 

Clearly there are uncertainties associated with truncation of the DGLAP 
evolution at NLO. Now, the DIS coefficient functions are known at NNLO [4]. 
Also valuable, almost complete, information has been obtained about the 
NNLO splitting functions [5]; indeed a range of compact analytic functions 
exist that are all compatible with this information [6]. We therefore study the 
effect of data cuts at NNLO. In going from NLO to NNLO the stable domain 
(xcut > 0.005, Ql ut > 7 GeV 2 , W^ ut > 15 GeV 2 ) has not increased as much as 
we might have hoped. Nevertheless there are advantages in going to NNLO [3]. 

To illustrate the value of having 'conservative' parton sets 1 at both NLO 
and NNLO we consider, as important examples, the determination of as(M^) 
from DIS data and the prediction of the cross section for W production, aw, at 

1 These sets are available at http://durpdg.dur.ac.uk/hepdata/mrs.html 
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the Tevatron. The values of as found in the NLO and NNLO global fits that 
produced the conservative sets of partons are given in Table 2 corresponding 
to the MRST03 entry, together with other recent determinations from DIS 
fits. The quoted errors reflect the tolerance A% 2 used in the various analyses. 
Remarkably, the determinations of as(M§) have converged approximately to a 
common value, even though they are based on different selections of the data. 



A X 2 



«s(Mf) 



± expt ± theory ± model 



NLO 

CTEQ6 
ZEUS 
MRST03 
HI 

Alckhin 



100 0.1165 ±0.0065 

50 0.1166 ±0.0049 ±0.0018 

5 0.1165 ± 0.002 ±0.003 

1 0.115 ±0.0017 ± 0.005 t °o 9 5 

1 0.1171 ±0.0015 ±0.0033 



NNLO 

MRST03 
Alckhin 



5 0.1153 ± 0.002 ±0.003 
1 0.1143 ±0.0014 ±0.0009 



Table 2: The values of a<?(M§) found in NLO and NNLO fits to DIS data. 
The experimental errors quoted correspond to an increase A\ 2 from the best 
fit value of x 2 . CTEQ6 [7] and MRST03 are global fits. HI [8] fit only a subset 
of the Fl p data, while Alekhin [9] also includes Ff and ZEUS [10] in addition 
include xF£ data. 



Fig. 2 shows values of the W production cross section (times the leptonic 
branching ratio Bi v = 0.1068) at the Tevatron energy y/s = 1.96 TeV, ob- 
tained from NLO and NNLO global fits of data, subject to various values of 
x cut and Q 2 ut . We see that the NNLO predictions are much more stable to 
variations of x cut - Note that at NNLO the conservative parton set predicts a 
small decrease of 0.7% relative to the NNLO default prediction (with x cut = 0, 
Q 2 ut = 2 GeV 2 ), while at NLO there is an increase. The conservative partons 
thus show greater convergence with increased perturbative order, than the de- 
fault predictions. The NNLO conservative value Bi„<j w = 2.67 nb, with an 
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Figure 2: Predictions for W production at the Tevatron, y/s = 1.96 TcV, for 
various values of Xcut , and for the conservative sets of partons (shown by the 
open symbols). 



expected total theoretical and experimental uncertainty of about 2%, may act 
as a very good luminosity monitor at the Tevatron. 
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